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Abstract 

We study the production of K^,p,u! and cp mesons in p C reactions on 
the basis of empirical spectral functions. The high momentum, high removal 
energy part of the spectral function is found to be negligible in all cases close to 
the absolute threshold. Furthermore, the two-step process {pN ttNN; ttN 
N + , p,uj, (f)) dominates the cross section at threshold energies in line with 
earlier calculations based on the folding model. 
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1 Introduction 



The production of heavy mesons in proton-nucleus reactions even below the free nucleon- 
nucleon threshold is of specific interest [|I], 0, ^, ||, |] as one hopes to learn about 
cooperative nuclear phenomena, high momentum components of the nuclear wavefunc- 
tion or in-medium modifications of the mesons themselves. Furthermore, a precise 
knowledge of the meson production channels in p + A collisions is a necessary step to- 
wards a microscopic understanding of meson production in nucleus-nucleus collisions, 
where in-medium modifications should be enhanced due to the higher baryon densities 
achieved j^, 0. Especially the production channels for the vector mesons is of par- 
ticular interest because their abundancy can be controlled by dilepton spectroscopy 
independently. A possible 'softening' of the p meson at high baryon density as 
possibly indicated by e~^e~ and data 0, |TU|, |TT|, ^ should also show up in the 

p° production on nuclei and could be controlled by the tt+tt" invariant spectrum. 

So far, the production of t] and mesons on nuclei has been addressed, both 
theoretically [T|, |14|, O, |l|] and experimentally [T|, H, |2l|, |2H|, H , whereas 



only a few model studies on the production of K and p,uj,4) mesons are available |^ 



|25[] . Most of the former approaches - though partly using the experimental momentum 
distribution for the target of interest - employ a quasi-free dispersion relation for the 
struck target nucleon and thus come in conflict with the A-body energy-momentum 
conservation close to the absolute threshold because the struck target nucleon is off- 
shell. 

Experimental information about the nucleon momentum and energy distribution 
is available from (e, e'p) reactions |^ . The data show in addition to the 'mean-field' 
spectral function a high momentum, high energy component which might be attributed 
to nucleon-nucleon correlations with high relative but low center-of-mass momenta [ p7|| . 
The influence of this 'correlation' component on heavy meson production in proton- 
nucleus collisions is of particular interest in our present study as well as the role of 
primary and secondary (pion induced) reaction channels. 

Our paper thus is organized as follows: In Section 2 we describe our approach for 
meson production in proton induced reactions with spectral functions and discuss its 
ingredients, i.e. the spectral function itself, the elementary production cross sections 
in pN and vrA^ collisions, the effective number of collisions as well as meson reabsorp- 
tion effects in Section 3. By discarding possible in-medium effects we then present in 
Section 4 our results for p C reactions for K~^, p, uo and mesons as a function of 
the beam energy Tq with their channel decompositions. A summary and discussion of 
open problems concludes this paper in Section 5. 



2 Reaction model for proton-nucleus collisions 



In terms of a multiple scattering theory the particle production in proton-nucleus col 
lisions can be described as an incoherent sum of two-body interactions 



28, 



30 



However, the individual interactions in nuclear matter might differ from those in free 
space due to the average interaction of the collision partners with the surrounding 
nuclear medium |l31 |. 
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The first term from tlie multiple scattering theory we denote as the 'direct' or 
'one-step' production process. Furthermore, neglecting the medium effects on the re- 
action partners, this approach corresponds to the 'impulse approximation' which can 
be applied whenever the relevant energy scale is large compared to the nuclear binding 
energy. Within the single-scattering approximation the Lorentz-invariant differential 
cross section for meson production in p + A collisions then can be written as 

d^PM ~ — T- 7/ ^ d^ 

with ^/s being the invariant energy of the incident proton and the struck nucleon 
defined as 

s = (E*+^f -(p*+q)2 (2) 

where Eq, Pq and u, q are the total energy and the momentum of the beam and target 
nucleon, respectively. Here the *-indices denote the influence of the average nuclear 
optical potential Uopt which changes the kinetic energy of the incident proton as 

Tq = To — Uopt, (3) 

where Tq is the beam kinetic energy and Eq = Tq + rriN, while m^v is the nucleon 
mass. Note that the momentum of the incoming proton pg also has to be modified 
properly in case of vector potentials. The momentum-dependent optical potential for 
cold nuclear matter at normal density po ~ 0.158/m~^ is about - 70 MeV for pq = 0, 
becomes repulsive at beam energies above 300 MeV and reaches an almost constant 



value of 30 - 50 MeV ^ ^ Thus the optical potential for the impinging 
proton is the first medium effect that has to be taken into account. 

The primed indices in Eq. (|l|) denote the meson momentum and energy in the cms of 
the beam and the individual target nucleon, while E'j^jd^apN^MX I d^p'j^j is the Lorentz 
invariant elementary cross section for meson production from p + N collisions in free 
space. The factor N^ff accounts for the effective primary collision number as well as for 
the in-medium absorption of the meson when propagating through the nuclear target 
as described in Section 3.4. In Eq. (|I|) the index Q stands for the Pauli blocking of the 
final states. 

The function Sp{q, relates the momentum of the target nucleon q with its energy 
uj. The free space relation, uj^ = q"^ + m^, is no longer valid for a system of bound 
nucleons because they are off-shell. In the mean-field approach the average total energy 
of the bound nucleon is [0, 

j2^'^j = ^-i^N- e) (4) 
i=i 

where Ma denotes the mass of the target while e stands for the binding energy per 
nucleon. Realistic spectral functions Sp{q,uj), as used in the present calculations, will 
be discussed in Section 3.1. 



-^We adopt the normalization / d^q dcj Sp{q,uj) — A, where A denotes the target mass. 
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As proposed in p^, pH| the dominant contribution to subthreshold heavy meson 



production should arise from a two-step reaction mechanism with an intermediate 
pion. In accord with |T3|, O] the cross sections for K^,p, uo and 0- mesons from the 



secondary pion induced reactions are calculated as 
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(5) 



where the double prime indices denote the system of the intermediate pion and a target 
nucleon, while the single prime indices are those for the beam and the first target 
nucleon. Moreover, ET,d^apN-,nX / d^Pn stands for the 7r-meson differential production 
cross section, which is calculated in line with (|I|), while atot is the total proton- nucleon 
cross section. We note that the on-shell elementary vr-meson production e.g. can be 
described within the isobar model The factor g^, in Eq. (j^) accounts for the 



probability that the pion interacts again with a target nucleon (cf. Section 3.4). 



3 Ingredients of the model 
3.1 The spectral function 

The in-medium relation between the energy uj and the momentum q of the target 
nucleon relevant to subthreshold particle production can be defined in analogy to Guet 



and Prakash |^T|, where the energy Ep^od (available for the meson production) is related 



to the incoming energy by total energy conservation as 

E* + Ma = Eprod + (m1„i + q2) + Eexc, (6) 

where Ma and Ma-i are the masses of the initial and residual nuclei, respectively, q is 
the internal momentum and i^exc stands for the excitation energy of the final nucleus. 
The excitation energy E^.^^. now is related to the removal energy Eji as [H2 



E,^, = Er + Ma- Ma-1 - rriN, (7) 

where Eji is the energy necessary to extract the nucleon from the nucleus. Thus 

UJ = rriN - Er (8) 

with Er > 0. The probability to find a nucleon with momentum q and removal energy 
Er in a nucleus now is given by the nucleon spectral function S{q,ER). Accordingly, 
the function Sp{q,uj) entering our calculations in (|l]),® can be evaluated from the 
spectral function S{q,ER) by substituting Eq. (H). 



The evaluation of the spectral function itself is quite an involved problem Q ^ 
and still a matter of discussion. In the present calculations we adopt S{q, Er) from 
Sick et al. [4^ obtained within the orthogonal correlated basis approach as well as a 
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phenomenological spectral function from Ciofi degli Atti and Simula |27|, that both 
provide a good description of (e, e'p) reactions. 

The S{q, E^) for from ref. is shown in Fig. |I] as a function of the nuclear 
nucleon momentum q and removal energy Er. In the upper part of Fig. 2 we display 
the spectral function up to Er k, 60 MeV on a larger scale in order to show the 'mean- 
field' distribution dominated by the Si/2 and P3/2 single-particle levels at Er ^ 35 MeV 
and Er ^ 20 MeV, respectively. As a 'background' under the horizontal shell structure 
one can notice a ridge roughly along 

Er ^ 45MeV - (9) 
2m N 

leading to a total available energy for particle production 

uj^rriN- AbMeV + (10) 

2'mN 

This is close to the dispersion relation often used in folding model calculations. 

The full spectral function up to Er = 500 MeV is given in the lower part where the 
'mean-field' contribution only shows up as tiny ridges, that extend up to g ~ 0.3 GeV/c, 
i.e. up to the region of the Fermi momentum. Beyond this value the distribution is 
dominated by the correlations. The dashed line indicates the maximum of the spectral 
function for high momenta attributed to short-range and tensor correlations as 

ER = 2e+-^ (11) 
2m]\i 

with e ~7 MeV The simple relation ( |TI| ) demonstrates that the high internal mo- 
menta correspond to high removal energies, too. This is just opposite to the behaviour 
for momenta below the Fermi level and will be quite important for subthreshold parti- 
cle production, because the contributions from the high momentum component of the 
nuclear wave function should be negligible close to the absolute threshold since the 
available energy (including Er) is below threshold. 

The relation between the high momentum and high removal energy was for the first 



time illustrated by Ciofi degli Atti, Pace and Salme considering the saturation of 
the momentum sum rule, i.e. 

^q)= / dERS{q,ER), (12) 

JEth 

where Eth is the single nucleon removal threshold. In Fig. ^we show the function ( [T^ ) 
for a ^^C target when integrating up to E^ax = 30, 100, 200 and 500 MeV. It be- 
comes clear that momenta below the Fermi momentum are essentially related to a 
small removal energy, which is again in reasonable agreement with the mean-field ap- 
proximation. 

The function (|12|), integrated up to infinity, is the nuclear momentum distribution. 
Respective momentum distributions for ^He,^^ C, and "^^^Pb, that have been calculated 
with the spectral functions from [P7|, are shown in Fig. R|. We display separately the 
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contributions from the uncorrelated part (dashed hne) as well as the saturation of 
the momentum sum rule for the correlated part at Emax='iO MeV (dotted line) and 
100 MeV (dashed-dotted line). The uncorrelated part of the spectral function includes 
the ground and one-hole states of the residual nucleus. The short-range and tensor 
parts of the realistic NN interactions deplete the states below the Fermi level and 
partially occupy the states above thus generating the component of S{q, E^i) with high 
momentum and high removal energy. The momentum distribution corresponding to 
the correlated part for Emax = 30 MeV (dotted lines) is found to be small for all nuclei 
in comparison to the 'mean-field' contribution. 



3.2 Cross sections for pion induced reactions 

For our calculations we will need the elementary particle production channels over 
a wide kinematical regime and, in particular, down to the reaction threshold, where 
the experimental data often do not exist. In this section we, therefore, formulate a 
parametrization that allows us to perform an extrapolation of the cross section into 
the unmeasured region. 

The cross sections for meson production from n + N collisions can be separated 
into two parts. The first one is an exclusive cross section related to binary processes, 
which are dominant at energies close to the reaction threshold, whereas the second part 
contains the inclusive production at higher energies. 

The cross section for the exclusive meson production from tt + N collisions has 
been calculated by Tsushima et al. E^; we adopt their results for our present study. 



The cross section for p, uj and 0-meson production from the reaction 

n + N ^ M + N (13) 

is most conveniently described within a resonance model. Assuming that the squared 
matrix element is proportional to a Breit-Wigner function the cross section for the 
reaction (|13]) can be written as 

.(.,N^ M ,N;s)^ ^g£±l f^f;^^^ X ,14) 

where J, Bin and B^ut are the resonance spin and the branching ratios of the incoming 
and outgoing channels respectively, while the factor R2 stands for the phase-space 
volume of the final particles, 

R2 = Ti \^l\s,mN,mM)/s^/\ (15) 

with tun and uim denoting the masses of the nucleon and meson, respectively. In 
Eq. (p^) k"^ is given by 

fc^ = A (s, rrtTT, mA?) (16) 

with the Kallen function 

^{.z,x,y)= z~{x + yf z - {x - yf /Az. (17) 
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The parameters and T from Eq. (|n|) are the mass of the 'baryon resonance' and full 
width, respectively. We note that within the resonance description we neglect coherent 
contributions from the available baryons as well as possible interference terms since we 
do not have experimental information on their actual magnitude. 

Within our semiphenomenological approach we now fit the experimental data on p, 
uj and 0-meson production in pion induced reactions in order to extract the mass and 
width of the 'effective' resonance R. Our fits to the data |^ are shown in Figs. 



with the parameters listed in Tab. |l|. For convenience, the data are plotted as a function 
of a/s — y/sth, with s denoting the squared invariant energy of the colliding particles 
and ^/sth = 2m N + niM in accordance with the reaction threshold. In Figs. §,|^J^ the 
fitted resonance cross section for the exclusive reactions Tt + N^uj + N^p + N^ip + N, 
are given by the dashed lines in comparison to the data from EH 



Since there are not enough experimental data available for the inclusive p, uj and 
0-meson production from it + N collisions that allow to construct reliable parametriza- 
tions, we calculate the cross sections for the inclusive meson production within the 
Lund-String-Model (LSM) from |^. The calculated results within the Lund-String- 



Model can be fitted by a function of the form 

ai^Tt+p-^ M + X)=a{x-lfx-^ (18) 
from threshold up to a few 100 GeV, where the scaling variable is defined as 

x = s/sth- (19) 

The parameters for all interesting channels a, b, c and Sth are listed in Tab. |^. 

Our fits to the LSM results are shown by the solid lines in Figs. ^,|,^and reasonably 
well reproduce the available inclusive data on vector meson production from inclusive 
processes at high energies. For practical purposes we will use the maximum of the 
'resonance' cross section (|T^ and the LSM parametrization (p!8| ) for the inclusive vector 
meson cross sections. 

The Lorentz invariant differential cross sections entering Eqs. (|I]) and (^, finally, 
due to two-body kinematics near threshold are completely determined when assuming 
an isotropic distribution in the pion-nucleon cms. 

3.3 Proton induced reactions. 

The cross sections for K^, p, u and 0-meson production from p + p collisions within 
the Lund-String- Model (LSM) are shown in Figs. as a function of the invariant 

collision energy ^/s in comparison to the experimental data from |^8|]. The solid lines 
show the fit to the LSM results according to Eq. (|TB]) with the parameters from Tab. ^ 
whereas in case of u production the full dots stand for the LSM results. We additionally 
show the cross sections for exclusive meson production within the One-Boson Exchange 
Model (OBEM) by the dashed lines in comparison to the respective data from |^ 
(triangles). Whereas in case of K~^, p, uj production the inclusive LSM results smoothly 
match with the OBEM calculations for the exclusive channels close to threshold, this is 



no longer the case for 0-meson production (cf. Fig. [T0|) because the string fragmentation 



model requires the formation of two ss-pairs, which shifts the threshold up to higher 
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energies. Since there are presently no experimental data available for y/s < 3.5 GeV 
in case oi p + p reactions we will adopt the maximum of the cross sections from the 
OBEM and LSM calculations for the inclusive 0-meson production as in the case of 
the + N reactions. 

Though our present results for the p + N channels have some uncertainty close to 
threshold, this will not show up sensitively in p + A reactions since the pion induced 
channels contribute more effectively (cf. Section 4). 



3.4 The effective collision number 



The factor Nf.ff{E) appearing in Eqs. (|T]), @ accounts for the A-dependence of the 
K~^, p, uo and 0-meson production as well as for the Final State Interaction (FSI) of the 



mesons in the nuclear medium. As first proposed by Margolis the A-dependence 
for particle production and propagation in the nuclear medium can be described by 
extending the Glauber multiple scattering theory |52|]. In accordance with | 



54 this 



gives 



N, 



eff 



+ 00 



bdb 



p{b, z)dz / d(j) 



2n 



/ rz r+oa 

expi^-apNj p{b,i)di- auN p(r[C])c/C 



where p(r) is the single-particle density in coordinate space (taken from 
to the target mass number. Here, y[Q is 

r[C] = ro(6,0,^) + Ce, 



(20) 
normalized 

(21) 



where b and z stand for the impact parameter and the 2;-component of the coordinate 
along the beam-axis, respectively, and e is a unit vector in coordinate space in the form 



e : = 



{sindcos(j), sin6sin(j), cos6), 



(22) 



while 6 is the meson emission angle in the laboratory system. Moreover, apN is the 
total p + N cross section while aMN is the M + N absorption cross section. 
Within the small angle approximation, i.e. for 6 ~ 0°, one can reduce (pO) to 



1 



eff 



O'pN — CTMN 



J d% {exp [-aMN TA{b)] - exp [-apN TAib)]) 



(23) 



with Ta(6) denoting the profile function of the target In Eq. ( p3D the second 

term is related to the effective number of nucleons (of the target) participating in the 
interaction with the impinging proton while the first term accounts for the final state 
reabsorption. 

As shown in ^ N^ff from Eq. ( p3D for aMN — can adequately be parame- 
terized by 

N,jf = A'-''^'-'\ (24) 

at bombarding energies above 0.8 GeV, where A is again the target mass number. 
Thus the A-dependence of the direct production mechanism is proportional to 
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whenever the produced meson only weakly interacts with the residual nucleus. In case 
of cTpTv ^ o'mn "we can - assuming a uniform distribution of the nuclear density - reduce 
the absorption term from (^) to the simple expression 

K = exp [-(TmnPoRa)] , (25) 

where po = 0.168 fm~'^ and Ra are the average density and the radius of the nucleus, 
respectively. We note that the factor k is frequently used to describe the reabsorption 
of the produced particles in finite nuclei in order to account for weak reabsorption 
effects. 

However, the application of ( [25| ) in case of p and u production is not adequate 
because the reabsorption cross section aMN appearing in (|20| - p5D is no longer small 
compared to (Jp^. In order to demonstrate the sensitivity of iVg// to the various ap- 
proximation schemes (^OD, (|23D , (^) we show in Fig.|lT] the effective collision number 



for ^^C as a function of the reabsorption cross section within the limits of Eq. (20) 
(solid line), Eq. (^) (dashed line) and Eq. ([25| ) (dotted line) which deviate substan- 
tially already for absorption cross sections of about 30 mb. Consequently, we will use 
Eq. (^) in the following to describe the inclusive production and absorption of mesons 
in p + A reactions with the absorption cross sections described in Section 3.5 

We, furthermore, note that by Eq. ( |20D we only account for reabsorption and ne- 
glect elastic meson-nucleon rescattering. This is of no special interest when considering 
inclusive (integrated) cross sections only as in Section 4, but will be essential when 
comparing the Lorentz invariant spectra with experimental data sets in a narrow kine- 
matic regime. For the latter purpose transport simulations as described in ref. 
will provide a more adequate description. 

The calculation of N^ff for the two-step reactions with an intermediate pion has 
been described by Vercellin et al. |^| within the Glauber formalism and is a more 



involved task. However, as shown in a practicable approximation is to use A^e// 



from Eq. (^) and to include the corrections from the two-step mechanism by a factor 
9n as ^ 

9. = 1 - + . (26) 

which describes the average probability for a pion to scatter again in the target with 
radius Ra- In case of small final meson {K^, p, u, 0) reabsorption the A-dependence 
of the inclusive production cross section due to the reactions with an intermediate pion 
roughly scales as 

N,ffXg^c:,A''^'-\ (27) 



3.5 



-nucleon cross sections 



The total cross sections for V+N interactions at high energies can be evaluated with the 
Vector-meson Dominance Model (VDM) from p, uj and 0-meson photoproduction 
or be obtained within the additive quark model relations [^, ^ from measurable cross 
sections. The results obtained with the VDM are in good agreement with those 
computed via quark sum rules, resulting in apN = ct^jn — 25 mb and a^N — 12 mb at 
vector meson momenta above 1 GeV/c. 
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The cross sections for p + N, 00 + N and (f) + N interactions in case of low relative 
momenta can be extracted from the vr + reactions by assuming that this is the 
dominant absorption channel. Detailed balance provides the relation 

a^+N^M+N _ {'^In + 1)(2/Af + 1) A(s, ruN, mu) .^g) 
aM+N->iT+N (2/^ + 1)(2/Ar + 1) A(s,mjv, m^) ' 

where s is the squared invariant energy, m^v, m^,/ are the masses, while Jat, and 
Im are the spins of nucleon, pion and vector meson, respectively. Here, the function A 
is given by (0). With our parametrizations for the cross sections Ut^j^n^m+n (0) 
then obtain the cross sections for vector meson-nucleon interactions via (ES). 



The cross sections extracted from detailed balance are only related to the particular 
inelastic reaction channel MN — > vrA^, and do not saturate the total cross sections even 
at momenta of 100 MeV/c. Indeed the 2TtN , SttN, rrniN or KKN final states as well 
as elastic scattering can also occur, which all contribute to the total cross sections. 
Since there are no reliable calculations on the + cross sections at low energy 
we adopt the results from the VDM. 



4 Results for p C 



After specifying all the elementary cross sections in proton and pion induced channels 
as well as the spectral function S'(g, p7| , ^ we now present the results for meson 
production cross sections in p C collisions according to (|I|) and (^). Though 
our primary interest are the vector mesons, we first discuss our calculations for 
production in comparison to the data from Koptev et al. |T^] and the results from the 
folding model [|13|. In Fig. |1^ we thus compare the cross sections within the spectral 
function approach for the one-step channel (dashed line) and the two-step mechanism 
(solid line) with the one-step (dash dotted line) and two-step result (dotted line) from 
the folding model with elementary cross sections as described in ref. [^ . 

Note that the folding model, contrary to the spectral function approach, is a pure 
mean-field approximation, where the nucleon momentum distribution for the target 
nucleus is calculated in the Hartree-Fock limit corresponding to the dashed line in 
Fig. ^. Thus the high momentum components, that are associated with the correlated 
part of the spectral function (cf. Fig. 3), are missing in the folding model calcula- 
tions 1 38]. The dispersion relation for a target nucleon in the folding model was taken 

uj{q) = 



as 



g2 + m% 



(29) 



where e is again the average binding energy per nucleon. In spite of the quite different 
assumptions in both approaches, the results for the i^'"'"-meson cross sections for one- 
step and two-step reaction channels are remarkably close to each other. 

Whereas at energies above 0.95 GeV both results are practically the same, the 
folding model with the quasi-free dispersion relation ( PUf ) for the target nucleon over- 
estimates the energy available in the collision and thus leads to slightly higher cross 
sections than our spectral function approach for both reaction steps at lower bombard- 
ing energy. Thus the underestimation of the data below about Tq = 0.92 GeV might 
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be considered as an evidence for additional many-body effects or further production 
channels that become important close to the reaction threshold. On the other hand, 
the dominance of the two-step reaction mechanism is the same within both approaches 
and can be considered as model independent. 

Since in the spectral function from we cannot separate the correlated and un- 



correlated parts, we adopt the spectral function from [27|, that allows the separation 
of the two components. Within the latter parametrization we can study the relative 
contribution from the correlated part of Sp{q, w) for meson production in p -\- A reac- 
tions. Our calculated results for i^"'"-meson production in p C collisions are shown 



in Fig. |I3| for the two-step mechanism (solid line) and primary reaction channel (dashed 
line) in comparison with the experimental data [|l9l. In line with our arguments pre- 
sented above the contributions from the correlated part of the spectral function indeed 
are small compared to the data. 

For more pedagogical reason we also discuss the results for the K~^-meson produc- 
tion cross section when employing the momentum distribution $(g), as obtained from 
integrating the spectral function from |^ over the removal energy, but using the free 
dispersion relation 



u = + q2, (30) 

for all momenta with m^r being the mass of free nucleon and q the nuclear momentum 
taken in accordance with the function $(g) from ||2^ shown in Fig. |^. 

In this unrealistic limit the one-step production mechanism (dashed-dotted line 
in Fig. 13) still underestimates the experimental data by more than a factor of 5, 
however, the results for the two-step mechanism now overestimates the data by almost 
an order of magnitude. This comparison clearly demonstrates that the high momentum 
component appearing in the spectral function approach cannot be exploited in a quasi- 
free production mechanism. 

We continue with our calculations for the vector mesons in p C reactions. The 



results are shown in Fig. [T^; here, the solid lines indicate the contribution from the 
two-step reactions with an intermediate pion (^, while the dashed lines correspond to 
the direct production mechanism (|I|) using again the spectral function from For 
all mesons (p, lj, 0) the two-step mechanism is clearly dominant close to threshold and 
kinematically is due to the effect that the Fermi motion of the nucleons is exploited 
twice in the two-step channels. The absolute cross sections at about 2 GeV are in 
the order of 10 /ib for p, cu-production, which is high enough to allow for a detailed 
experimental study at COSY or CELSIUS. 



5 Summary 

In this work we have presented an approach for meson production in p + A reactions on 
the basis of direct and two-step production channels (including an intermediate pion) 
on the basis of the nuclear spectral function without incorporating any selfenergies 
for the mesons. The 'elementary' cross sections for p, uj and <^ production have been 
evaluated within the LUND string formation and fragmentation model and within 



a resonance model for vrA^ collisions close to threshold, respectively. The available 
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experimental data from [|8| for the exclusive and inclusive reactions are described 
sufficiently well, however, our extrapolations for the various channels close to threshold 
should be controlled experimentally in the near future. 

Whereas in previous reaction models only the target momentum distribution has 
been employed using a quasi-free dispersion relation for the target nucleons, our present 
formulation is based on the full spectral function S{q, uj) where the nucleon may be far 
off-shell and high momentum, high removal energy components appear as well. 

This approach was first used by Debowski, Grosse and Senger ^ for proton 
induced particle production. While those authors applied the method to their own 
measured differential spectra for i^^-meson production on and ^°®P6, we have 
tested our approach for K'^ production in comparison to the experimental data from 
Koptev et al. [|19| using the spectral function for ^^C from ref. . 



What emerges from a detailed study of the spectral function is that one should use 
a quasifree dispersion relation (|T0|) only for nuclear momenta below the Fermi level. 
For higher values of g a dispersion relation with (^TJ) should be used for the origin of 
the high momentum component accounting for two-body correlations. 

It is clearly seen that the two-step mechanism with an intermediate pion gives by 
far the dominant contribution at subthreshold energies as in the folding model 
but the actual cross section for i^"^-meson is underestimated below about 920 MeV 
indicating the need for further reaction channels or many-body effects. This result 
contrasts with the folding model where a better description of the data is achieved. 
We attribute this to the fact that the folding model - adopting a quasi-free dispersion 
relation for the target nucleon - overestimates the available energy in nucleon-nucleon 
collisions when approaching the absolute reaction threshold. 

On the other hand, when integrating the spectral function over the energy En, 
we obtain a momentum distribution showing very high momentum components that 
can be attributed to short-range nucleon-nucleon or tensor correlations. However, these 
high momentum components on average are associated with high removal energies, too, 
such that they cannot be exploited for meson production close to the absolute reaction 
threshold. Thus, meson production close to threshold does not provide information on 
the high momentum components in finite nuclei or short-range two-body correlations. 

Within our approach we have also calculated the cross sections for the vector mesons 
p, CO, (p in p C collisions from threshold up to a few GeV. We obtain estimates for 
their production in nuclei without incorporating any selfenergies for the mesons. These 
cross sections provide a baseline for experiments at COSY or CELSIUS that address 
the in-medium properties of the vector mesons as well as their final state interaction 
with nuclei. Here, in context with chiral symmetry restoration, the p and tu-mass 
at normal nuclear density po is expected to be reduced by about 18% @], which will 
sizeably enhance their production cross section as compared to our calculations for free 
vector mesons presented above. A detailed study of their production cross section as a 
function of Tq and target mass A might help to separate the effects from the real and 
imaginary part of the meson selfenergies. 
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Table 1: The effective mass and width of the baryonic resonance in Eq. (|T^, as well 

as B = Bin X Bout X (2 J + 1). 



Meson 


Mn (GeV) 


r (GeV) 


B ifihGeV-'') 




1.809 


0.99 


413 


CO 


1.809 


0.99 


302 





1.8 


0.99 


5.88 



Table 2: The parameters of Eq. (|T8D for pion induced reactions. Here the index 
imphes a shift x —>■ x — 1.3. 



Meson 


sth (GeV2) 


a (mb) 


h 


c 


P" 


2.917 


3.6 


1.47 


1.25 




2.958 


4.8 


1.47 


1.26 




3.831 


0.09 


2.54 


2.1 



Table 3: The parameters of Eq. ( |T8|) for proton induced reactions. Here the index '^^^ 
implies a shift x ^ x — 1.3. 



Meson 


Sth (GeV2) 


a (mb) 


h 


c 


K+ 


6.49 


1.12 


1.47 


1.22 


P° 


7.01 


2.2 


1.47 


1.1 




7.06 


2.5 


1.47 


1.11 




8.38 


0.09 


2.54 


2.09 


7] 


5.88 


2.5 


1.47 


1.25 
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Figure 1: The spectral function S{q,Eji) for ^^C from [H6|. The dashed hne shows the 



relation (pAl). The contour lines decrease by a factor of 1.5 from line to line. 

Figure 2: Saturation of the momentum sum rule for ^^C. The lines indicate the func- 
tion (0) integrated up to energy E^ax=^0 MeV (dashed), 100 MeV (dotted), 200 MeV 
(dash-dotted) and 500 MeV (solid line). The arrow indicates the Fermi momentum. 



Figure 3: The momentum distribution for ^He, ^^C, and "^^^Pb according to p^. The 
dashed lines result from the uncorrelated part of the spectral function. The dotted 
lines indicate (|1^) for the correlated spectral function integrated up to Emax = 30 MeV 
and the dashed-dotted lines up to Emax = 100 MeV. The solid lines display the sum of 
the uncorrelated and correlated parts for E^ax oo. 

Figure 4: Experimental cross sections for the exclusive reaction tt^ + p + p 

(triangles) and the inclusive reaction tt^ + p ^ + X (full dots) from [HH]. The lines 
represent our parametrizations as discussed in the text. 

Figure 5: Experimental cross sections for the exclusive reactions tt"*" + n uj + p 
(squares) and 7i~ +p oj + n (triangles) and the inclusive reaction +p — > uj+X (full 
dots). Experimental data are taken from The lines represent our parametrizations 
as discussed in the text. 

Figure 6: Experimental cross sections for the exclusive reactions 7r+ + n (f) + p 
(squares) and it~ + p —>■ (p + n (triangles) and inclusive reaction it~^ + p —>■ + X (dots) 
from |H8[. The lines represent our parametrizations as discussed in the text. 



Figure 7: Cross sections for the exclusive reaction p + p ^ + A + p (triangles) 



and the inclusive reaction p + p — > + X (dots) from . The solid line shows the 



calculation within the LSM, while the dashed line refers to the OBEM. 

Figure 8: Cross sections for the exclusive reaction p + p p° + p+p (triangles) and the 



inclusive reaction p + p p° + X (dots) from |48]. The solid line shows the calculation 



within the LSM, while the dashed line refers to the OBEM. 



Figure 9: Experimental cross sections for the exclusive reaction p + p uo + p + p 
(triangles) from |^ in comparison to the OBEM calculation (dashed line). The full 



dots represent the results from LSM for the inclusive reaction p + p ^ uo + X, while 
the solid line shows our fit with Eq. 16. 

Figure 10: Cross sections for the exclusive reaction p+p — >■ (triangles) and the 



inclusive reaction p + p ^ (j) + X (dots) from . The solid line shows the calculation 
within the LSM and the dashed line our results within the OBEM for the exclusive 
channel. 
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Figure 11: The effective collision number calculated with (pOD-solid, (|23D -dashed and 
with (^)-dotted line as a function of the meson absorption cross section a in mb. 



Figure 12: Cross section of i^'"'"-meson production for C as a function of the beam 
energy Tq. The solid and dashed lines show the calculations with the total spectral 
function for the two-step and direct reaction mechanism, respectively, employing the 
spectral function from The dotted and dashed-dotted lines show the calculations 
with mean-field approach (dashed curve in Fig. 3) for the two-step and direct reaction 
mechanism, respectively, employing a quasifree dispersion relation. The dots show the 
experimental data from |]rU|. The solid and the dotted lines were calculated with the 
same elementary cross sections. 



Figure 13: Cross section of fT^-meson production for C as a function of the beam 
energy Tq. The solid and dashed lines are results using only the correlated part of the 
spectral function from |^ for the two-step and direct reaction mechanism, respectively. 
The dotted and dashed-dotted lines show the on-shell calculations for the two-step and 



direct reaction mechanism, respectively. Experimental data are from |19 



Figure 14: Cross section for p , uo and 0-meson production in p C reactions as 
a function of the beam energy Tq. The solid and dashed lines show the calculations 
with the total spectral function |46] for the two-step and direct reaction mechanism, 
respectively. 
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